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ABSTRACT 

We present a complete atlas of the Cygnus Loop supernova remnant in the 
light of [0 III] (A5007), Ha, and [S II] (AA6717, 6731). We include low-resolution 
(25") global maps and smaller fields at 6" resolution from observations using 
the Prime Focus Corrector on the 0.8-m telescope at McDonald Observatory. 
Despite its shell-like appearance, the Cygnus Loop is not a current example of a 
Sedov- Taylor blast wave. Rather, the optical emission traces interactions of the 
supernova blast wave with clumps of gas. The surrounding interstellar medium 
forms the walls of a cavity through which the blast wave now propagates, 
including a nearly complete shell in which non-radiative filaments are detected. 
We identify non-radiative shocks around half the perimeter of the Cygnus 
Loop, and they trace a circle of radius R = 1?4 (19 pc) in the spherical cavity 
walls. The Cygnus Loop blast wave is not breaking out of a dense cloud, but is 
instead running into confining walls. Modification of the shock velocity and gas 
temperature due to interaction of the blast wave with the surrounding medium 
introduces errors in estimates of the age of this supernova remnant. The 
optical emission of radiative shocks arises only where the blast wave encounters 
inhomogeneities in the ambient medium; it is not a consequence of gradual 
evolution to a global radiative phase. Distance measurements that rely on this 
uniform blast wave evolution are uncertain, but the radiative shocks can be used 
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as distance indicators because of the spherical symmetry of the surrounding 
medium. The interstellar medium dominates not only the appearance of the 
Cygnus Loop but also the continued evolution of the blast wave. If this is a 
typical example of a supernova remnant, then global models of the interstellar 
medium must account for such significant blast wave deceleration. 



Subject headings: ISMdndividual (Cygnus Loop)-shock waves-supernova 
remnants 



1. Introduction 

Supernova remnants greatly determine the large-scale structure of the interstellar 
medium. The energy of supernova remnants heats and ionizes the interstellar medium 
(ISM), and their blast waves govern mass exchange between various phases of the ISM. In 
doing so, supernova remnants (SNRs) influence subsequent star formation and the recycling 
of heavy elements in galaxies. Global models of the interstellar medium that include a 
hot ionized component ( |Cox fc Smith 1974j ; [McKee fc Ostriker 1977| ) are sensitive to the 



supernova rate, the persistence of their remnants, and the sizes they attain. A simple 
calculation of the last of these assumes that the blast wave expands adiabatically in a 
uniform medium once the blast wave has swept up mass comparable to the mass of the 
ejecta. During this Sedov- Taylor phase, the radius of the SNR as a function of E51, the 
initial energy in units of 10 51 erg, n , the ambient number density in units of cm -3 , and £4, 
time in units of 10 4 yr, is R = ^(E^x/no) 1 ^ 5 ^ 5 pc in a medium where the mean mass per 
particle is 2.0 x 10~ 24 g. This phase will last until radiative losses become important. The 
beginning of this subsequent phase, marked by the initial loss of pressure behind the blast 
wave, occurs at t — 1.9 x 10 A El( 14: n~ A ^ 7 yr, when the radius is R = 16.2El( 7 n~ 3 / 7 pc flShull 
fc Draine 1987| ), although the radiating shell is not fully formed yet. 



We approach these large-scale questions with analysis of complete images of a particular 
supernova remnant, the Cygnus Loop, in three optical emission lines. This supernova 
remnant appears to be a limb-brightened shell at radio (|Keen et al. 1973|) , infrared ([Braun 



fc Strom 1986|) , optical flFesen, Blair, fc Kirshner 1982Q , and X-ray (|Ku et al. 198g ; 



Levenson et al. 1997] ) energies, which at first glance suggests that it is presently in the 



transition to the radiative stage. The Cygnus Loop has the advantages of being nearby, 
bright, and relatively unobscured by dust. This allows us to examine in detail the evolution 
of various portions of the shock front and to determine physical parameters, such as shock 
velocity and local ambient density, as they vary throughout the remnant. Despite its 
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appearance, the Cygnus Loop is not a current example of blast wave propagation in a 
uniform medium at any stage. Instead, its evolution is governed by the inhomogeneous 
interstellar medium, which we map using the shock as a probe. 

Many of the features we discuss have been noted by others. Oort (1946) first suggested 
that the Cygnus Loop is an expanding supernova shell. Spectroscopy of radiative shocks 
in selected locations (e.g., Miller 1974j , |Raymond et al. 1980a| , and Fesen et al. 1982) 
combined with theoretical models of these shocks (e.g., |Cox 197% popita 1977| , [Raymond 



1979| , and (Shull fc McKee 1979| ) has been used to derive the physical conditions of the 



observed shocks. We utilize the radial velocity measures of Minkowski (1958), Kirshner & 
Taylor (1976), Greidanus & Strom (1991), and Shull & Hippelein (1991) to discern some 
of the three-dimensional structure that is ambiguous from the data we present. Many 
non-radiative or Balmer-dominated shocks in the Cygnus Loop have been identified (e.g., 
Kirshner fc Taylor 1976 . Raymond et al. 1980b , [Trcffcrs 1981| , Fesen et al. 1982, Fesen 



Kwitter, fc Downes 1992] , and [Hanson fc Winkler 1992j ). Our observations qualitatively 



match these, and we rely on these works and others ( [Raymond et al. 1983| ; |Long et~aT| 
1992; Hester, Raymond, &: Blair 1994) for quantitative measures of parameters such as 



shock velocity and preshock density. McCray & Snow (1979) and Charles, Kahn, & McKee 
(1985) have suggested that the Cygnus Loop is the result of a cavity explosion, and we 
adapt this global model to interpret the surrounding interstellar medium, as well. 

This paper is a companion to the soft X-ray survey presently in progress with the 
ROS AT High Resolution Imager ( Graham fc Aschenbach 1996| ; [Levenson et al. 1997|) . With 



these two surveys, we examine the Cygnus Loop as a whole, not restricting our investigation 
only to those regions that are exceptionally bright or that appear to be particularly 
interesting. We hope to understand both the global processes and the specific variations 
that are responsible for the emission we detect. The X-rays probe hot (temperature T ~ 10 6 
K) gas that shocks with velocities v s ~ 400 km s" 1 heat. The optical emission is expected 
from slower shocks (v s ^ 200 km s -1 ) in which the post-shock region cools to temperatures 
T ~ 10 4 K, yet the most prominent regions at optical wavelengths are also bright in X-rays. 
McKee & Cowie (1975) suggested that the broad correlation of X-ray and optical emission 
is the result of a blast wave propagating in an inhomogeneous medium. In this scenario, 
the shock is significantly decelerated in dense clumps of gas, while portions of it proceed 
unimpeded through the lower-density intercloud medium. We apply the principles of this 
basic cloud-blast-wave interaction to a range of locations in the Cygnus Loop. In particular, 
we refine the cavity model introduced in Levenson et al. (1997), using these optical data to 
constrain the current ISM in the vicinity of the Cygnus Loop and to determine how the 
stellar progenitor modified it in the past. 
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We present the observations in §2. We describe them in detail, noting individual 
regions of interest, and we use these data to measure the physical conditions of the blast 
wave and the ambient medium in particular locations in §3. The purpose of the detailed 
examination is to combine the results in a complete map of the surrounding ISM. We 
present this three-dimensional model while providing a coherent explanation of the history 
that accounts for it in §4. We predict the future of this SNR and relate its fate to more 
general theories of supernova modification of the interstellar medium in §5 and summarize 
our conclusions in 56. 



2. Observations 

We obtained narrow-band images of the Cygnus Loop in 1995 August and September 
and 1995 September on the 0.8-m telescope at McDonald Observatory with the Prime Focus 
Corrector (PFC; |Claver 1992| ). The PFC system is a 5-element catadioptric corrector located 
at the F/3 prime focus of the telescope and projects a flat focal plane onto a Loral- Fairchild 
2048 x 2048, 15/im-pixel charge-coupled device. The camera and corrector system has a 0'.'3 
spot size in the V-band and projects each pixel to l'/35 on the sky, producing a 46' x 46' 
field of view. The interference filters each have a rectangular transmission profile and full 
width at half maximum FWHM = 40A in the converging beam of the PFC. The filter 
centers in the F/3 beam are 5007A, 6723A, and 6563A to detect line emission from [O III] 
(A5007), [S II] (AA6717,6731), and Ha+[N II], which we refer to throughout as Ha. Each 
individual integration of approximately 600 s duration was bias-subtracted, flat-fielded, and 
scaled by its median sky value before being combined with other images. Every region of 
the Cygnus Loop was observed at least twice through each filter in order to remove cosmic 
rays and exclude bad regions of the detector, and most regions were observed three times 
in Ha. We obtained a total of 62 observations in [O III], 89 in Ha, and 54 and [S II]. 

Over the large spatial scale of the Cygnus Loop (3° x 3?5), each observation is a 
distinct projection of the sky onto the plane of the detector, so images must be remapped 
to a common projection before being co-added. An astrometric solution for each integration 
was calculated based on the Hubble Space Telescope Guide Star Catalog. Alignment does 
not rely on overlapping observations of a particular region but depends instead on the 
astrometry of individual integrations. The final images are the result of pixel-by-pixel 
averaging in the combined image plane using nearest-neighbor matching and rejecting bad 
regions of the detector. The global maps of the Cygnus Loop with 25" pixels are shown with 
linear scaling in Ha, [O III], and [S II], and in a false-color combination with logarithmic 
scaling (Figure []]). While the observations at different wavelengths resemble each other 
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broadly, their variation is apparent in the combined image. In this display, blends of Ha 
and [S II] alone are yellow, [O III] and [S II] together appear as magenta, and Ha and [O III] 
combine to make cyan. Only the regions that are displayed as white exhibit strong emission 
at all three wavelengths. In addition to the global maps, we present 5 smaller fields, each 
95' across, at 6" resolution (Figs. @-|6|). Although we first describe the current physical 
conditions in specific regions of interest, our ultimate goal is to combine the understanding 
of the particular sites into a complete account of the surrounding ISM, its modification by 
the progenitor star, and the passage of the blast wave through it. 



3. Emission Morphology and Structure of the ISM 

3.1. Observed Morphologies 

The radiative shocks that are responsible for most of the optical emission from the 
Cygnus Loop result in two characteristic morphologies: filaments and diffuse emission. 
Figure [7] identifies examples of these features and other structures to which we will refer 
by name. The morphological distinction is the result of different viewing geometry. As the 
"wavy sheet" model (Hester et al. 1983; Hester 1987| ) describes, bright, sharp filaments are 



the result of long lines of sight through tangencies of the shock front. When the shock front 
is viewed face-on, the emission appears more diffuse. 

Combined with the X-ray data, these observations provide more specific information 
on the evolutionary state of the Cygnus Loop and the distribution of inhomogeneities 
of the ISM around it. The widespread X-ray emission indicates that the global shock 
velocity is several hundred kms -1 , so the optical emission comes from regions of denser gas 
that have decelerated the blast wave. The correlation of the X-rays with optically-bright 
regions make clear that this SNR is not globally in the gradual transition to the radiative 
phase that is the result of blast wave propagation in a uniform medium. The blast 
wave has recently and suddenly decelerated where radiative shocks are detected. This 
rapid, environmentally-imposed evolution is distinct from the simple model prediction for 
incremental change in a homogeneous medium. 

These complete data sets demonstrate that many characteristics that have been 
investigated in small regions of the Cygnus Loop are widespread, and we adapt the 
physical explanations of these limited regions in previous work to understand the Cygnus 
Loop globally. For example, the bright emission from radiative shocks is found to be 
stratified at select observed locations ([Fesen, Blair, fc Kirshner 1982|) , and we show this 
stratification to be typical of all prominent radiative shocks. The [O III] extends toward 
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the projected exterior, and the Ha and [S II] are located toward the projected interior. 
The [S II] is farthest to the interior and is particularly diffuse. The stratification is a 
function of the column density that the shock has swept up. The incomplete [O III] 
emission is concentrated behind shocks that have encountered relatively little gas, while 
where the shock lags because it has swept up a greater column of gas, Ha and then [S II] 
are prominent. The strong [O III] emission requires that the shock has progressed through 
a column density N ~ 3 x 10 17 cm" 2 , while the complete recombination zone develops 
after the swept-up column N ^ 10 18 cm -2 for v s < 150 km s ([Raymond et al. 1988|) . 
The apparent thickness of these zones and the extensive regions of bright emission in a 
combination of all the observed lines is a consequence of projection effects. The curved 
surface of the blast wave encounters the cloud at different projected radii, so some lines of 
sight intersect multiple stages of shock evolution. 

Around most of the perimeter of the Cygnus Loop, filaments are detected in Ha alone. 
These Balmer-dominated filaments are the result of non-radiative shocks ( |Chevalier fc 



Raymond 1978| ; |Chevalier, Kirshner, fc Raymond 1980Q . Unlike the more common radiative 



shocks, the Ha emission here does not come from the extended cooling zone behind the 
shock. Instead, it is the result of collisional ionization before excitation in the immediate 
post-shock region, or the downward cascade toward the ground state following charge 
exchange. In either case, neutral atoms are required in the post-shock region. Because these 
do not survive long, the emission originates at the shock front. Thus, we use the Balmer 
filaments to trace the shock front directly. This Ha emission is intrinsically weak, so it is 
most obvious when the line of sight through the shock front is long. This is the case at the 
projected edge of the SNR, where the line of sight is tangent to the spherical blast wave. 



3.2. Projected Edges 

The projected edges of the Cygnus Loop are the simplest regions to understand because 
the geometric complications are minimized there. The observed structures may be slightly — 
but not significantly — in the foreground or background and still appear near the projected 
edge. The prominent regions, NGC 6992 at the northeast (a = 20 h 56 m 24 s , 5 = 31°43'; 
Fig. 0) and NGC 6960 at the west (a = 20 h 45 m 42 s , 5 = 30°43'; Fig. |), are examples of 
interactions of the blast wave with large clouds ( [Hester fc Cox 1986| ; Hester et al. 1994; 



Graham et al. 1995| ; |Levenson et al. 1996| ), and the corresponding X-ray enhancement 



arises from reflected shocks that propagate through the hot, shocked SNR interior. 



The observed length scales of optical and X-ray emission from NGC 6992 and NGC 
6960 reveal the sizes of these large clouds. At the 770 pc distance of the Cygnus Loop, the 



- 7- 



coherent networks of filaments over 20' imply that the clouds have lengths of about 10 19 
cm. Balmer filaments, high-resolution X-ray observations or both constrain the location of 
clouds along the line of sight. X-ray emission extends beyond both of these bright optical 
edges, where some portion of the three-dimensional blast wave is not yet impeded by the 
clouds. Thus, the clouds do not extend over the complete edge of the blast wave along the 
line of sight. 

In contrast, the southeast knot at a = 20 h 56 m 20 s , S = 30°25' (Fesen et al. 1992; 
Graham et al. 1995|) is an example of a cloud that is clearly extended along the line of 
sight but is not necessarily large across the plane of the sky. The nearly-circular edge of the 
X-ray SNR and the Balmer filaments are distinctly concave at the southeast knot (Fig. |]a), 
and there is no X-ray emission to the exterior. Thus, this cloud must be at least 11 pc long 
along the line of sight in order to impede the entire projected edge of the blast wave at this 
location. Again, the presence of the complete optical cooling structure indicates that the 
shock has progressed through a column N ^ 10 18 cm~ 2 here. 

The breakout to the south, centered on a = 20 h 50 m ,5 = 29°10' and extending 
approximately 1° across (Fig. ||), is the most significant departure from circularity in the 
remnant. The greater advance of the blast wave in the breakout compared with the circular 
loop north of S = 29°40' indicates that the ambient density is lower in this southern region. 
Balmer-dominated filaments surround much of the breakout, where the blast wave has 
now encountered atomic gas. There are sections of incomplete shocks to the east and west 
(e.g., a = 20 h 53 m ,5 = 29°20' and a = 20 h 48 m , 5 = 29°15') and fully radiative shocks at 
the east (a = 20 h 52 m ,5 = 29°15'). Emission is absent across the face of this area. In the 
plane of the sky, the breakout is not very extended. Either the blast wave has only recently 
entered this low-density region, so it has not had enough time to advance significantly, or 
the extent of the low-density material is limited, so the blast wave no longer continues to 
progress rapidly. The stratification of emission around the edge of the breakout is apparent 
on smaller scales than those observed elsewhere in the Cygnus Loop, and the breakout 
exhibits smaller regions of common overlap. In the northeast, for comparison, the distinct 
regions of Ha and [O III] are typically 3' wide, and they are observed together over 10' 
scales, whereas in the breakout, the exclusive and common zones are both typically 1' wide. 
These characteristics are consequences of the more limited scale of interactions with clouds 
in the breakout. There is no evidence for indentation or blast wave progress beyond the 
optical-emitting regions, which could be due to either intrinsically small clouds or early 
stages of encounter. In either case, at least some dense material surrounds the SNR at its 
southern extremity. 

The exterior of the southwest quadrant exhibits many filaments that are nearly aligned 
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with one another (Fig. [5]). Some of these are part of the radiative shock structure related to 
the southern end of NGC 6960, while farther away from the center, these are non-radiative 
and incomplete shocks (Fesen et al. 1992). Along some of the filaments, the transition 
from non-radiative to incomplete shock structure is clear where [O III] emission arises. The 
multiple filaments are easily understood as tangencies of a wavy sheet. These filaments 
are shorter than those at the north and east of the Cygnus Loop. This indicates that the 
ISM is clumpier here, which would distort the blast wave in many places and account for 
the number of filaments that we observe. Shock diffraction through low-density (Sp/p ^ 1) 
regions may be responsible for the large curvature of some of these filaments, as Fesen et 
al. (1992) suggest, although multiple interfering shock fronts are not necessary to produce 
the different filaments. They arise from multiple projections of the same shock front. 



3.2.1. Balmer Filaments and X-ray Emission 

The Balmer-dominated filaments around the periphery of the Cygnus Loop indicate 
that the fraction of atomic gas surrounding the SNR is high. The extreme circularity of 
these filaments requires that the blast wave remain nearly spherical with little deformation 
at the projected edge. With the exclusion of the breakout region, we fit them with a circle 
of radius R = 1?4, which is equivalent to 19 pc at a distance of 770 pc. In Figure the 
model circle is drawn on the soft X-ray image obtained with the ROSAT High Resolution 
Imager ( Levenson et al. 1997|) . We calculate the correlation of radial variation as a function 
of angular separation and the corresponding Fourier coefficients for the Balmer filaments. 
Excluding the breakout, the amplitude of the m = 2 mode is less than 0.02, which implies 
that the surroundings of the Cygnus Loop are homogeneous on large scales. 

The individual filaments tend to be very long, typically extending over more than 
40'. The continuity and smoothness of the Balmer filaments indicate the uniformity of 
the medium in which these shocks propagate. For constant ram pressure pv 2 s behind 
the shock, variations in density and velocity are related Sp/p = —25v s /v s . In terms 
of the radius, R, and time, t, v = rjR/t, where rj is of order unity; for free expansion 
rj = 1, and for an adiabatic supernova remnant in a uniform medium, rj — 2/5. Thus, 
5p/p = (—2/r])(5R/ R)(t/ 5t) . The observed radial variations of a particular filament 
are SR/R ~ 0.01. An approximate timescale of these shocks is 1000 yr, and the age of 
the Cygnus Loop is about 14,000 yr (cf. §4.2]), so the overall density variations in the 
present medium of the Balmer-dominated filaments are of order 5p/p ~ 0.3r]. We expect 
2/5 < r) < 1, so 0.1 £ 5p/p £ 0.3. 

The Balmer-dominated filaments define the current location of the blast wave and 
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mark the presence of neutral material. Detailed studies of particular locations at the 
northeast have used these non-radiative shocks as density probes ( [Raymond et al. 1983 



Long et al. 1992|; Hester et al. 1994) and derive densities n ~ 1 cm 3 . The total observed 



Balmer line intensity is sensitive to density, shock velocity, and geometric perspective. All 
of the filaments around the periphery share the same edge-on geometry. Because of the 
remnant's circularity, it is reasonable to assume that the shock velocity is nearly constant 
around the extreme projected edge except at the breakout. The filaments that have been 
studied in detail to determine shock velocity and ambient density are the brightest ones, 
however, so these regions are most likely in the transition to incomplete radiative emission, 
and these densities are expected to be higher than elsewhere. Thus, we conclude that where 
typical non-radiative filaments are observed around the periphery, the density is constrained 
n cm -3 , and most likely n ~ 0.1 cm -3 . 

We roughly calculate the Ha surface brightness of a non-radiative shock in the Cygnus 
Loop. We examine a portion of the blast wave in the southeast, near a = 20 h 56 m , 5 = 30°00'. 
This is not an exceptionally bright filament, and there is no obvious interaction with a large 
cloud here. The whole filament extends over more than 20', so the blast wave edge is smooth 
here. We estimate a surface brightness of 4.2 x 10 -6 erg s -1 cm -2 sr -1 from this filament. 
Assuming that this line of sight passes through two folds of the spherical blast wave and the 
ISM does not distort it, the intensity of a single shock surface is 2.1 x 10 -6 ergs -1 cm -2 sr -1 . 
The line intensity, /, depends on the energy of the transition, hi/, and is given by 

hv q ex _ 1 _ 2 —1 
1 = — nnVg — ergs cm sr , 
4tt qi 

where nu is the preshock neutral density and the ratio of excitation and ionization rates 
Qex/Qi ~ 0.2 for Ha ([Raymond 199 1| ) . If we reasonably constrain the shock velocity 
200 < v s < 400 km s -1 , these observations limit the density 0.3 > nn > 0.15 cm -3 at the 
location of the non-radiative shock. The lower velocity limit is based on observations of 
Balmer filaments ( [Raymond et al. 1983 ; |Long et al. 1992 ; Hester et al. 1994), which are 



likely to represent the slowest non-radiative shocks, as noted above, and the upper bound 
assumes that the forward shock does not heat the X-ray-emitting medium to temperatures 
greater than 2 x 10 6 K. 

The Balmer filaments are related to adjacent X-ray emission: the optical filaments 
form the exterior boundary of low surface brightness, limb-brightened X-rays QLevenson et 



al. 1997| ). This characteristic morphology reveals where the blast wave is decelerated in 



the cavity walls. We combine the density determined from the Balmer observations with 
X-ray temperature measurements to calculate the unperturbed blast wave velocity and 
the original density of the cavity, and we compare this with the cavity density based on 
the observed X-ray surface brightness of the interior of the SNR. A bimodal temperature 
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distribution has been detected in the Cygnus Loop, with peaks at 1.5 x 10 6 and 5 x 10 6 
K in EXOSAT observations with the channel multiplier array (sensitive over 0.05-2 keV) 
and the medium energy experiment (sensitive over 1-50 keV; [Ballet fc Rothenflug 1989| ). 
We suggest that the higher temperature component is due to the reflected shock and that 
previous low single-temperature measures are contaminated by the slow, forward shock in 
the cavity walls. We model these simply as a reflected shock and ignore non-equilibrium 
effects, which are probably important. The physical parameters temperature, pressure, and 
density, are related by the equations of conservation of mass, momentum, and energy flux 
across shock fronts. Assuming the temperature of 1.5 x 10 6 K corresponds to the decelerated 
shock in the cavity wall and T = 5 x 10 6 K corresponds to the reflected shock, the current 
blast wave velocity is around 300kms _1 , and the original, unperturbed blast wave velocity 
was around 500 km s -1 . 

The reflected shock model also predicts that the original cavity density, n c , is related 
to the shell density n s /n c = 5. Combined with the surface brightness measurement of the 
Balmer filaments, this yields an atomic shell density n s = 0.2 cm -3 and a cavity density 
n c = 0.04 cm -3 if the shell is entirely neutral. As described below ( §4.1| ), the shell is 
partially ionized, so the original cavity density was n c = 0.08 cm -3 . In addition to neglect 
of non-equilibrium effects, other sources of error include confusion of the reflected shock 
and the unperturbed shock in the cavity, uncertainty of the derived X-ray temperature, 
and uncertainty of the Balmer intensity of a single shock surface derived from observation 
through the spherical edge of the blast wave. 

The predicted density is consistent with that derived from X-ray observations of the 
projected interior. In a featureless region of the southeast near a = 20 h 55 m , 5 = 30°12', 
the median HRI surface brightness is 0.01 counts s _1 arcmin -2 . We compare this with a 
model of a Raymond-Smith plasma at T e = 3.5 x 10 6 , correcting for absorbing column of 
Nh = 6 x 10 20 cm~ 2 . Assuming the emitting region is 22 pc deep (the predicted line of 
sight through the SNR at this radius), the observed count rate corresponds to a density 
n e = 0.2 cm -3 . This is the density of the shocked material; the original cavity density would 
have been n c = 0.05 cm -3 . 



3.2.2. Photoionization 

One unusual feature at the east is a region of very smooth emission adjacent to the 
bright structure of NGC 6995 at a = 20 h 57 m 30 s , 5 = 31°00'. This area is obvious in Ha and 
detectable in [S II] but absent in [O III]. There are no X-rays here ( Levenson et al. 1997| ), so 



a fast shock wave has not passed through. The X-ray map shows a clear indentation in this 
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region, where dense gas has significantly impeded the blast wave over large scales, similar 
to the southeast cloud. There is no other direct evidence, such as a Balmer filament, that 
the blast wave is presently located toward the exterior of this region. We propose that this 
gas is photoionized by a shock precursor, and it has not been shocked. 

Radiation from the shock-heated gas of the interior of the Cygnus Loop serves as 
the ionization source. Typically, when the shock velocity v s ^ 100 km s _1 , much of 
the subsequent radiation is at UV wavelengths ( [Raymond 1979| ; |Shull fc McKee 1979| ; 



Popita et al. 1984]) , which can ionize the surrounding gas, both ahead of and behind the 
shock front. The hardness of the ionizing spectrum increases with shock velocity, and 
the resulting emission from photoionized gas is a function of shock velocity and preshock 
density ( [Sutherland, Bicknell, fc Dopita 19§5| ; popita fc Sutherland 1996] ; [Morse, Raymond,] 



|fe Wilson 1996] ). The absence of [O III] in the ionized precursor indicates that the shock has 
intermediate velocity (v s ^ 200 km s -1 ), which is consistent with UV spectroscopy of Blair 
et al. (1991) who find v s ~ 170 km s -1 is typical in this region. The bright Ha emission 
ahead of the shock front is probably a consequence of relatively high density; the Balmer 
luminosity of the precursor scales linearly with density in the steady case (Raymond 1979 



jShull fc McKee 1979[ Popita fc Sutherland 1996| ). A second important factor is that the 



shock is fully radiative, having been decelerated in the dense cloud, so the shock front does 
not move through the ionized region before the ions have time to radiate (cf. [Raymond et 



al. 1988|). Finally, the column density ahead of the shock must be great for the precursor 



emission to be observable. 

The images do allow us to estimate crudely the density of the photoionized material. 
We assume that the length of the photoionized region along the line of sight is equal to 
its extent north-south in the plane of the sky and use the observed surface brightness to 
predict a density n ~ 60 cm~ 3 . Spectral data would allow a more exact determination of the 
characteristics of the emitting region and quantitative measurements of the photoionizing 
source. 

The blast wave history causes the ambient medium to be photoionized here and 
not elsewhere. The penetrating X-rays of very fast shocks (y s ^ 400 km s -1 ) do 
not efficiently ionize the surrounding medium; the optical depth to 200 eV photons 
T200eV = N/(l.l x 10 20 cm~ 2 ). Only the softer UV-dominated spectra of slower shocks are 
effective ionizing sources. In other regions of the Cygnus Loop, the recent shock velocity has 
been too great to produce much photoionization, and there has not been enough time for 
the ionizing field of the post-shock cooling and recombination zone to develop. Thus, most 
of the surrounding medium has remained sufficiently neutral to exhibit Balmer-dominated 
filaments. The exceptional photoionization at the eastern edge is consistent with the cavity 
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model of this SNR. The blast wave has not slowed sufficiently elsewhere for most of the 
shocked gas to be a useful source of ionizing photons. Over nearly all of the Cygnus Loop, 
the deceleration has been recent, as the blast wave has run into the walls of the cavity. 
Only in this region has the softer spectrum of a slower shock existed long enough to ionize 
the preshock medium so that its radiative transitions can be detected. 



Photoionized regions have been observed around other SNRs, including N132D ( [Morse 



Winkler, fc Kirshner 1995| ) and N49 QVancura et al. 1992| ) in the Large Magellanic Cloud, 
for example. One difference between the Cygnus Loop and the former SNR, however, is 
that while in the Cygnus Loop the photoionized precursor is the exception, detected only 
in a limited region, in N132D it is the rule, occurring over the majority of the younger 
remnant's perimeter. Another difference is that the exceptional photoionization ahead 
of the Cygnus Loop blast wave is adjacent to the region that contains its brightest soft 
X-ray knots ([Levenson et al. 1997 ), while the photoionization and X-ray surface brightness 
variations are not correlated in N132D QMorse et al. 19961) - These differences are likely 
to be a consequence of the different sources of ionizing photons and the density of the 
surrounding medium. First, in N132D, the X-ray-emitting shock can account for only half 
of the ionization ( [Morse et al. 199~6|) , so another significant source must also contribute to 
the ionization in this case. The ubiquitous radiative shocks in clouds and O-rich filaments 
that Morse et al. (1996) propose as this additional source are not preferentially associated 
with the X-ray emission. Second, a large molecular cloud extends over the entire southern 
hemisphere of this SNR ( [Banas et al. 19971) . While the present blast wave is not in the 
dense core that is observed in CO, the outer regions of the cloud are expected to be 
over-dense with respect to the typical ambient ISM. The higher densities that make the 
ionizing precursor visible are therefore characteristic of the majority of the perimeter of 
N132D. In contrast, we suggest that it is the unusual conditions of extreme blast wave 
deceleration and related high density that result in both the bright X-ray emission and the 
photoionized precursor in this small area at the eastern edge of the Cygnus Loop. These 
conditions make the Cygnus Loop precursor similar to that of N49. This middle-aged SNR 
is associated with a small molecular cloud on its eastern limb flBanas et al. 199"7[ ). The blast 
wave interaction with the cloud results in local optical emission and X-ray enhancement, 
blast wave deceleration observed in a range of shock velocities ( |Vancura et al. 1992j ), and a 
prominent photoionized precursor on the eastern limb of N49. 
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3.3. Projected Interior 

The interior of the Cygnus Loop contains diffuse emission and filaments. Both 
morphologies share some characteristics, such as exclusive [O III] emission (e.g., 
a = 20 h 48 m ,5 = 30°30' in diffuse regions and a = 20 h 52 m ,5 = 30°30' in filaments), 
and blends of Ha and [S II] (e.g., a = 20 h 50 m ,5 = 29°45' and a = 20 h 50 m ,5 = 29°45'). 
A combination of emission in all three lines is observed only in diffuse emission 
(e.g., a = 20 h 52 m , 5 = 31°10'), and Ha alone is observed only in filaments (e.g., 
a = 20 h 50 m ,<5 = 30°30'). 

Projection effects are more severe in the interior, but these optical images and 
high-resolution spectroscopy ( [Kirshner fc Taylor 1976| ; |Greidanus fc Strom 1991| |Shull fe| 



Hippelein 1991|) together help unravel the tangled Cygnus Loop. In all cases, the optical 
emission is the result of the interaction of the three-dimensional blast wave surface with 
dense structures in the ISM. The emitting regions are truly at the outside of the SNR, 
but some appear projected to the interior. The brighter filaments tend to have lower 
absolute radial velocities and narrower velocity dispersions than the fainter diffuse emission 



( [Kirshner fc Taylor 1976| ; |Shull fc Hippelein 1991| ). This is consistent with the model in 



which filaments are long lines of sight through the edge of a wavy shock front, and diffuse 
emission results when the blast wave is observed nearly face-on. The shock velocity is 
primarily perpendicular to the line of sight in the former case and along the line of sight in 
the latter case. In some instances, adjacent filamentary and diffuse emission are related, 
being tangential and normal views of the same portion of the blast wave. Elsewhere, 
spatially separated shock regions are projected onto the same location in the observer's 
view, causing significant confusion. 

The extensive diffuse emission at the western interior near a = 20 h 47 m — 49 m and 
5 = 30°40' — 31°20' suggests that a large wall of dense gas is present across the entire face of 
the SNR (Figs. |3| and The irregularity of the diffuse emission is due to the dumpiness 
of the cloud. This emission is blueshifted with respect to nearby gas QGreidanus fc Strom| 
|1991| ; |Shull fc Hippelein 1991|) so it is on the front face. The bright filaments of the "carrot" 



(at a = 20 h 49 m , 5 = 31°40'), having positive radial velocities, are on the rear face. 

The filaments at a = 20 h 49 m , 5 = 31°00' (Fig. |3|) are distinct on small scales, similar 
to the edges of the breakout, although here the region of common overlap is absent. The 
origin of the stratification is the same: the limited physical scale of the interaction with the 
intervening cloud, either because the cloud is intrinsically small or because the encounter 
is recent and the blast wave has not progressed far into the cloud. Along a particular line 
of sight, only a single coherent section of the blast wave is detected. Unlike the areas 
where strong emission is observed at all wavelengths, here there are not many projected 
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overlapping regions of various stages of shock evolution. The filaments that are bright in 
Ha and [S II] illustrate an earlier section of the cloud shock than do the [O III] filaments. 
In this area on the western half of the Cygnus Loop, there is no planar blast wave that 
actually moves to the east, however. Rather, the interaction with a cloud is progressing 
from west to east as the blast wave surface expands radially. The intervening cloud must 
be oriented along the line of sight to some degree, as shown in Figure |9|. 

Another region of diffuse emission and fainter filaments runs north-south across the 
center of the remnant, around a = 20 h 50 m — 52 m ,5 = 30°20'— 50' (Figs. | and |). The 
high-resolution spectra in the vicinity are of the diffuse component alone, which includes 
very bright central knots. These clumps are blueshifted, so these clouds are on the near 
face. The westernmost of these filaments, at a = 20 h 50 m 20 s , is detected only in Ha. 
The blast wave has run into an extended wall of gas, enabling us to look through a long 
edge of the shock front. These filaments are broader than most of the other pure Balmer 
filaments, which is expected in this scenario. The broad filaments appear where the blast 
wave is viewed through a point of inflection, whereas the sharper, narrow ones are viewed 
at a tangency. For these non-radiative filaments to be observable, the wall of gas must be 
extremely smooth at the present projected location of the shock front. The wall of gas is 
likely to be related to the source of the rest of the optical emission at the center of the 
SNR, not the clouds that are illuminated as the "carrot," although this is uncertain. 



4. History of the Cygnus Loop 

4.1. Modelling the ISM 

Except in the single photoionized section, the optical emission traces dense shocked 
regions. Using some knowledge of shock evolution and geometry, we map the surrounding 
ISM. The basic principles applicable to these observations are: (1) optical emission due 
to radiative shocks requires that the blast wave has swept up a significant column density 
(N ^ 10 18 cm~ 2 ), which is likely in regions of increased density (n ^ lcm" 3 ); (2) sharp 
filaments are edge-on views of the blast wave surface, and diffuse emission is the result 
of face-on views of this surface; (3) the optical line emission from a radiative shock is 
progressively dominated by [O III], Ha, then [S II]; and (4) filaments of Ha alone are the 
result of non-radiative shocks in material that has a significant neutral fraction. 

With these data, we construct a three-dimensional model of the interstellar medium 
around the Cygnus Loop . Over much of its surface, the blast wave is propagating through 
large clouds. Long filaments, detected in Ha around the perimeter and in other lines 
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across the interior, indicate that a spherical shell of atomic gas surrounds the SNR. The 
dense clumps of gas and the atomic shell together form the walls of a cavity in which the 
supernova occurred. Slices through the Cygnus Loop parallel to the plane of the sky from 



the near side to the far side (Fig. W) show some of the three-dimensional structure of the 
ISM before the supernova event. The center of the Cygnus Loop is at latitude b = —8°. At 
the rear of the cavity, toward the Galactic plane, is a large wall of gas that is detected in 



neutral hydrogen ( PeNoyer 1975|) and in the IRAS sky survey ( [Braun fc Strom 1986Q . The 



cloud that causes the bright northeastern emission is also on the far side. Clumps at the 
edge of the breakout and at the northwest are on the near side. Some clouds, including the 
southeast knot and the molecular gas at the west, extend along the line of sight and appear 
in all slices. 

The Cygnus Loop is not the result of a supernova explosion in an arbitrary region of 
the interstellar medium. Rather, the blast wave propagates through a medium that the 
progenitor star has processed. Several characteristics imply a massive progenitor, the sort 
most capable of modifying the surrounding ISM, which Charles et al. (1985) and Levenson 
et al. (1997) discuss. One strong piece of evidence is the paradoxical morphology of the 
Cygnus Loop. At optical and X-ray wavelengths, it appears to be nearly circular, yet it 
exhibits tremendous structure on smaller scales. The optical data alone could be interpreted 
as a SNR in transition from the adiabatic to the radiative phase, but the presence and 
correlation of X-rays require that overall the blast wave still has (or until very recently 
had) a high velocity, v s ^ 400 km s _1 . A model of blast wave interaction with large clouds 



accounts for this correlation (Hester fc Cox 1986; Graham et al. 1995; Levenson et al 



|1996|) , but such large clouds cannot be typical of the pre-supernova interstellar medium. 
If they were, the SNR would not appear so circular, and X-ray emission due to clouds 
evaporating in the interior would be observed. Thus, although large clouds are found around 
the periphery of the Cygnus Loop, they were not located within the current blast wave 
prior to the supernova event. This structure is expected in the vicinity of an O or B star. 
During its main-sequence lifetime, such a massive star homogenizes the surrounding ISM, 
creating a uniform region at constant density, temperature, and pressure. The star's UV 
emission removes clumps of denser gas within the H II region, either destroying the clouds 
through photoevaporation or relocating them to the exterior by the rocket effect ( [McKeeJ 



Van Buren, fc Lazareff 1984|) . The result is a homogeneous spherical region surrounded by 
dense clouds. A star of type earlier than BO will clear a cavity with a radius of about 50 
pc (McKee et al. 1984 ). The Cygnus Loop cavity is smaller, so its progenitor must be a 
later- type star flCharles, Kahn, fc McKee 1985| ), yet one still able to create the photoionized 
cavity. Thus, we suggest that a BO, M w 15Mq star was the progenitor of the Cygnus 
Loop. 
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The observer's line of sight through shocked clouds at the projected edge is preferentially 
longer than through those that appear at the projected center. Photoevaporation of clouds 
progresses radially outward from the central star, so lines of sight tangent to the cavity are 
biased to be through the longer dimensions of surviving clumps. Thus, optical emission on 
the projected interior is not as bright as on the projected edge of the SNR, despite their 
common cause. The same effect is even more pronounced in the X-ray data. 

The atomic shell in which the non-radiative filaments appear is also a consequence of 
the progenitor's evolution. While material is eroded from clouds during the star's main 
sequence lifetime, increasing the density of the H II region, the radius of the ionized region 
decreases as the outer portions recombine, creating the shell ( [Elmegreen 1976p . The atomic 
gas of this recombined region together with the cloud remnants on the periphery form the 
walls of the cavity (Fig. |TUp. The Cygnus Loop blast wave has expanded through the 



uniform cavity and is now hitting the surrounding walls QLevenson et al. 19971) . 



The atomic shell is not actually a sharp density discontinuity. It is more likely a 
gradual transition from the cavity interior to the maximum shell density. At the end of the 
progenitor's main sequence lifetime, the pressure in the former H II region declines as the 
gas cools and recombines. (See |Shull et al. 1985| for a discussion of this phenomenon applied 



to N49.) In the extreme case of a sudden loss of pressure at the boundary, a rarefaction 
wave would move back through the dense shell resulting in a density profile of the form 
n oc (1 + ^^) 3 flZeldovich fc Raizer 1966|) , where the distance, r, from the original shell 



edge is constrained Satnc < r < in terms of sound speed, a, and red giant lifetime, tnc- 
The drop in pressure is not so sudden, however. Over t^c ~ 10 6 for M ~ 15Mq ( |Maeder~%; 



Meynet 1989| ), about half the atoms in the former H II region have recombined. Thus, the 



density profile of the shell edge is likely to be intermediate between a sharp discontinuity 
and the offset power law above. 

The Cygnus Loop has been suggested previously to be the result of an explosion in an 
incomplete cavity. Both Falle Sz Garlick (1982) and Shull & Hippelein (1991), for example, 
argue for a partial cavity at a density discontinuity, such as the edge of a molecular cloud. 
In these models, only half of the blast wave runs into the cavity boundary. Consequently, 
the bright radiative emission of NGC 6992 and NGC 6960 and the corresponding X-ray 
emission are physically separated by great distances along the line of sight: the former are 
the result of radiative shocks in the cavity wall, and the opposite hemisphere of the blast 
wave propagates through the low-density medium outside the cloud to cause the latter. In 
contrast, the model presented here and in Levenson et al. (1997) has a completely defined 
cavity, the atomic shell and dense clouds together forming its spherical boundary. The blast 
wave does not break out of the edge of the molecular cloud into lower-density material; 
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instead, it slows down when it hits a wall of denser gas. Although NGC 6992 and NGC 
6960 and the X-rays that appear to the exterior of these radiative regions are two slightly 
different portions of the blast wave, they are not significantly separated. Furthermore, 
some of the large clouds are obviously extended along the line of sight, not restricted to 
the rear boundary (Fig. fLO)) . These optical and previous X-ray data QLevenson et al. 



|1997|) demonstrate that the cavity wall is complete. All around the projected edges of the 
Cygnus Loop, the Balmer filaments and X-ray enhancements mark the smooth component 
of the cavity wall. Even the region called the breakout illustrates where the blast wave runs 
into a boundary of dense gas, not continued expansion in a low-density medium. Only the 
eastern third of the projected interior of the Cygnus Loop lacks direct evidence for a cavity 
wall, although the smooth atomic shell that is apparent in Balmer-dominated filaments 
elsewhere would not be detected across the face of the SNR. Only clumps of gas could be 
observed, and their distribution depends on the arrangement that existed before the ISM 
was processed by the progenitor. The southeast of the Cygnus Loop is farthest from the 
galactic plane, so we expect the original filling factor of clouds to have been lowest there. 



Other SNRs, including RCW 86 QVink, Kaastra, fc Bleeker 1997] ) and N132D ( Morse et 



al. 1996| ), have been identified as cavity explosions. [Shull (1985) identifies N49 as a cavity 
supernova, but Vancura et al. (1992) argue that the dense ambient material associated 
with the single observed CO cloud accounts for the data better.] In these cases, interaction 
of the blast wave with the cavity wall enhances the observable properties of the SNRs. 
RCW 86 shares many features of the Cygnus Loop, including similar size, a limb-brightened 
(though very incomplete) X-ray shell (Vink et al. 1997), and Balmer-dominated filaments 
to the exterior of the X-ray emission ( (Smith 1997|) . N132D is a younger remnant, and the 
apparent cavity is smaller. The progenitor's wind during its main sequence lifetime may 
form the cavity boundary, or the scale of photoevaporation may have been limited in this 
dense environment adjacent to a molecular cloud ( [Banas et al. 199~7|) . The bright interior 
X-ray emission (|Morse et al. 1996]) does imply that the cavity interior was not uniform. 



4.2. Age and Distance 

Although the shell-like visual morphology of the Cygnus Loop and its X-ray 
limb-brightening broadly suggest that this is an example of a supernova remnant in the 
Sedov- Taylor or adiabatic expansion phase, in detail this model is inconsistent with the 
data presented here and in Levenson et al. (1997). The observed X-ray and optical 
limb-brightening are a consequence of interaction with the inhomogeneous interstellar 
medium, not the result of evolution of the blast wave in a uniform medium. Assuming 
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that the blast wave has slowed sufficiently during its progress through a uniform medium 
to produce optical line emission (i.e. to v s <^ 200 km s -1 ) yields an age estimate of about 
50,000 yr ( [Minkowski 1958] ; |Falle fc Garlick 1982] ). This is much too great because the age 



is inversely proportional to the blast wave velocity, and only in propagating through dense 
clouds has the shock been decelerated to v s ^ 200kms _1 . The Balmer filaments mark 
the location of the unimpeded, spherical blast wave, and the X-ray emission immediately 
behind these filaments demonstrates that the typical shock velocity is v s ^ 400 km s -1 , 
which is still to great for radiative cooling of the Sedov- Taylor shell to be significant. This 
SNR has not made the global transition to the radiative evolutionary phase. The X-ray 
emission alone has been used to determine an age of around 18,000 yr ( [Rappaport et af 



1974fc [Ku et al. 1984] ), assuming that this is presently an adiabatic blast wave expanding in 



a uniform medium; t = 2R/5v s , where R = 19 pc and v s = 400 km s 

The X-ray data provide a more accurate measurement of the age, but it is still 
likely to be an upper limit because these calculations contain uncertainties in the shock 
temperature and the velocity that is derived from it. Only the current unperturbed shock 
velocity is relevant to the age in this model, but any observation of the SNR includes 
contributions from a variety of temperatures, heated by different stages of the blast wave's 
history. Observing the projected edge of the SNR minimizes confusion due to overlapping 
projections of different temperature components, but the interaction of the blast wave with 
the inhomogeneous ISM greatly alters the shock velocity and the corresponding postshock 
gas temperature. 

When the blast wave hits a density discontinuity, as it does over most of the surface 
of the Cygnus Loop, the continuing shock is suddenly decelerated and a reflected shock 
propagates back toward the interior, further heating the dense, hot, shocked material. At 
a density contrast of 10, for example, the twice-shocked matter will have a temperature 
approximately 1.5 times the singly-shocked material at the same radius, which leads to an 
overestimate of the unperturbed shock velocity by about 20%. Both the reflected shock and 
the forward shock through the dense cloud enhance the X-ray surface brightness and thus 
bias the mean measured temperature to favor blast-wave-cloud interactions rather than the 
undisturbed blast wave. 



If the X-ray measurement of two temperature components ( [Ballet fc Rothenflug 



1989] ) due to the decelerated cloud shock and the reflected shock is correct and the 
original velocity of the unperturbed blast wave was closer to 500 km s" 1 (§§3), the 
corresponding age t = 14, 000 yr. A reliable determination of the age of the Cygnus Loop 
is required to determine the initial energy, E Q , accurately (§[[]). Adopting the parameters 
n Q = 0.08 cm"" 3 , R = 19 pc, and t = 14,000 yr, and a partially-ionized cavity interior, we 
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find E D = 2 x 10 50 erg. 

Because the radiative emission shows where the blast wave has run into clumpy material 
and decelerated, care is required to use it to measure the distance to the Cygnus Loop. 
Minkowski (1958) determined the expansion velocity to be 116 km s -1 , fitting an ellipse to 
radial velocity measurements as a function of radius. This, combined with Hubble's (1937) 
observed proper motion of filaments of 0.03" year -1 yields the frequently-cited distance 
of 770 pc. This calculation assumes that all the optical filaments are part of a common 
expanding surface. In fact, the optical emission reveals different regions of inhomogeneity 
in the ISM. A second complication is the distinct geometry of the diffuse and filamentary 
emission on the face of the SNR. One is the result of a face-on shock, and the other is an 
example of a shock viewed edge-on. The latter does not have the radial velocity structure 
as a function of projected radius that is assumed in the model calculation. A third source of 
confusion arises because the observed filaments are the cooling zone behind the shock front, 
so their apparent motion may be the result of further development of the postshock region, 
not motion of the same clump of gas. The cooling zone may be propagating, without any 
physical motion of a parcel of gas in the ISM. 

Despite these difficulties, however, the radiative shocked gas can be used to measure 
distance to the Cygnus Loop or other supernova remnants when the distribution of clumps 
in the surrounding medium is spherically symmetric. In this case, although the emission 
comes from distinct regions of the ISM, not a single expanding blast wave surface, the 
different clouds acquire the same velocity when the blast wave passes through them. This 
approach is not appropriate where the surrounding medium is asymmetric, however. IC 
443, for example, is asymmetrical, and its corresponding distance determinations are highly 
uncertain ( preen 1984| ). 

The non-radiative filaments can also be used to determine distance, as Hester, 
Danielson, & Raymond (1986) have done, finding a distance of about 700 pc. These 
structures are advantageous because proper motion and spatial velocity can be calculated 
from the same filament. There is some uncertainty in the shock velocity because the derived 
value depends on the amount of ion-electron temperature equilibration (Chevalier et al. 
1980). Measuring proper motion of filaments also introduces some uncertainty. All filaments 
are detected where the observer's line of sight through the shock front is long, which may be 
the tangent surface of a spherical bubble. If it is instead the deformed surface of the blast 
wave that appears as a filament, however, the apparent filament will not necessarily persist 
or be detectable over the timescales required to obtain proper motion measurements. Long 
baselines over which to measure proper motion are desirable, but the oldest observations 
are sensitive to the brightest filaments. These bright filaments are exactly the locations 
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where the shock is likely slowed in denser gas, so the present velocity differs from the spatial 
velocity that corresponds to the proper motion. With these caveats in mind, we adopt a 
distance of 770 pc to the Cygnus Loop. 

5. Consequences 

The Cygnus Loop blast wave is now in the walls of the cavity, the boundary between 
the low-density former H II region and the ambient ISM. Blast wave propagation in the 
atomic shell will be a long-lived phenomenon. The smooth material of the shell recombined 
as clumps were photoevaporated in the interior. We crudely calculate the thickness of the 
shell assuming a constant number of ions in the H II region in both the initial Stromgren 
sphere, in which only the less-dense intercloud material was ionized, and in the final 
Stromgren sphere of radius Rf, in which all material, having mean density n m , is ionized. 
The typical, unprocessed ISM determines the initial conditions; we assume the intercloud 
density rii C = 0.1 cm -3 , and a volume filling factor of 0.03 of clouds of density n c \ = 1 cm -3 , 
so n m = 0.13 cm~ 3 . Conserving ion number and neglecting increased recombination as 
clouds are evaporated, we obtain for the shell thickness 

R sh = R f ((n m /n iC ) 1/3 - 1), 

so using the observed radius Rf = 19 pc yields R s h = 1.5 pc. A blast wave of velocity 
v s = 300kms _1 will propagate through the shell for 5000 years. 

Eventually the blast wave will propagate through the denser medium of the cavity 
walls and clouds everywhere. It will rapidly decelerate, not only in the select locations 
where it has already encountered dense clumps that protrude into the cavity. Some of the 
surrounding material consists of large dense clumps like the ones observed presently in the 
very extended regions of radiative shocks. These, too, will light up and will be detectable 
at optical wavelengths. The southeast cloud may be an example of the earliest stages of 
interaction with a large cloud, and it may later resemble the northeast region flGraham 
et al. 1995|) . We also expect the entire northern limb, which in the optical regime is now 
detected primarily in bright Balmer filaments, to turn into an extended radiative region. 
The incomplete shocks along the northern limb are the early stages of interaction with a 
larger structure that is perhaps directly associated with the bright northeastern region. The 
entire western face of the Cygnus Loop will soon be bright when the blast wave has fully 
encountered the wall of gas that we have noted there. If the long filaments across the center 
of the SNR are the result of another wall of gas, the middle will also light up with optical 
emission in several hundred years. 
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Large clouds of gas will significantly affect blast wave evolution. The shock will no 
longer be adiabatic when the cooling timescale becomes comparable with the dynamical 
timescale. In large clouds, where density n c i = lcm -3 , radius R c i = 5pc, and in which the 
shock velocity v c i = 200 km s -1 , the dynamical time R c i/v c i ~ 2 x 10 4 yr. In these clouds, 
the cooling time 3kT ' j 'n^A(T) ^4x 10 4 yr, in terms of the cooling function, A, having units 
ergcm 3 s _1 . If the covering fraction of clumps in the surrounding medium is sufficiently 
great, they will force the Cygnus Loop into the radiative phase, and energy losses will 
further decelerate the blast wave. The time and radius at which this occurs are distinct 
from self-similar calculations, such as those of McKee & Ostriker (1977) and Shull & Draine 
(1987) in which the gradual transition to the radiative phase happens as the blast wave 
propagates through a uniform medium. 

We use the pressure-driven snowplow model ( [McKee fc Ostriker 1977| ) to set an upper 
limit on the final size of the SNR. In this model, work of the blast wave pushing into 
the ambient material and radiative losses of the hot interior are taken into account. In 
the Cygnus Loop, the radiative losses of the outer shell also will be significant, slowing 
the blast wave and making the final radius smaller. During this phase, PV 1 is constant, 
where P is the pressure inside the SNR, V is its volume, and 7 is the ratio of specific 
heats. For 7 = 5/3, this phase ends when R = Rpds(Ppds/ Pism) 1 ^ 5 , in terms of the 
external pressure of the interstellar medium, Pism, and the radius and pressure inside the 
supernova remnant at the onset of the pressure- driven snowplow stage, Rpos an d Ppos, 
respectively. At this final time, the pressure inside the supernova remnant is equal to the 
ambient pressure. For the Cygnus Loop, we take Pism = 10~ 12 dyne cm -2 , Rpns = 20 pc, 
and Ppds = 4 x 10~ 9 dyne cm -2 , which places an upper limit on the final radius R < 105 pc. 

Assuming that blast waves propagate through uniform media over great distances yields 
misleading conclusions about the size of SNRs and how they distribute energy through the 
ISM. Global models of the interstellar medium that include the hot contribution of SNRs 
( |Uox fc Smith 19"73| ; |McKee fc Ostriker 1977| ) are sensitive to the supernova rate, the final 



sizes of supernova remnants, and their persistence as hot bubbles. These factors determine 
whether the dominant component of the ISM is hot and in which smaller regions of warm 
and cold gas are embedded, or whether the warm and cold components form the substrate 
in which the distinct hot bubbles of SNRs are located. For the former to occur, the onset 
of the radiative phase of the SNR must be postponed until R > 100 pc, independent of 
its prior evolution ( |Cox 1986|) . (This calculation assumes that the supernova rate is 1 



per 30 yr in a disk of radius 15 kpc, yet the conclusion holds for a range of reasonable 
supernova rates.) Although the cavity of the Cygnus Loop is small because its progenitor 
is one of the least massive stars to undergo core collapse, the initial mass function favors 
these lower mass (and smaller cavity) progenitors. Furthermore, even more massive stars, 
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where cavities of 50 pc radius are expected, will be forced into the radiative phase and not 
expand to overlap if the density of the cavity walls is sufficient. The cavity walls are indeed 
likely to be denser on average than the typical ISM because the walls contain much of the 
material that was cleared out of the cavity interior. The conditions of the ISM in general, 
however, constrain the mass that will be commonly available to form clouds and cavity 
walls around pre-supernova stars. If the hot, low-density medium is pervasive, cavity walls 
typically would not be sufficiently massive to induce the radiative phase early, and SNRs 
would then normally attain larger radii before their blast waves decelerate. In any case, 
the typical environments of supernova progenitors must be considered. If dense, molecular 
material normally surrounds these stars, the radiative phase will be induced early, even if 
globally the ISM has low density. 

The cavity structure is significant also in determination of the mass of gas that a 
supernova remnant heats. The onset of shell formation for a supernova blast wave expanding 
in a uniform medium occurs at R oc n~ 3 / 7 , so the shocked mass is then M oc n~ 2 /"' ' . If the 
ambient density varies by three orders of magnitude, the heated mass changes by only a 
factor of 7, for example. In contrast, when the supernova environment is a cavity with 
dense walls, the initial size of the cavity constrains the mass of the hot interior. The cavity 
interior is rarefied, so the heated cavity mass is significantly lower. In the Cygnus Loop, 
assuming a cavity density n c = 0.08 cm" 3 , the heated mass is 20Mq, whereas in a uniform 
medium of this density, the heated mass would be greater than 400Mq. 

6. Conclusions 

This optical emission line atlas of the Cygnus Loop supernova remnant provides the 
information to render a portrait of the surrounding ISM. The ambient medium, which the 
massive stellar progenitor shaped as it evolved, consists of many large, dense regions. These 
will significantly affect the subsequent development of the blast wave. The blast wave will 
decelerate, and it will no longer have sufficient velocity to excite the material through which 
it propagates to X-ray-emitting temperatures. 

Although this is a study of a particular object, it has far-reaching consequences for the 
ISM as a whole. Supernova remnants provide the energy to heat the ISM, affect the velocity 
dispersion of interstellar clouds, and set the stage for future generations of star formation. 
Thus, any global model of the gas in the Galaxy critically depends on the evolution of SNR 
blast waves in their inhomogeneous environments, of which this example is typical. While 
optical emission-line characteristics identify supernova remnants, they may preferentially 
select those whose evolution the extant ISM determines. Other SNRs need to be examined 



-23- 



over a broad range of energies in the same careful way to draw more certain conclusions 
about their net effect on the interstellar medium. 

This research was supported in part by the Packard Foundation. 
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Fig. 1. — Complete narrow-band maps of the Cygnus Loop obtained using the PFC at 
McDonald Observatory. Each individual integration covered a 46' x 46' field using l'/35 
pixels. The mapping procedure used to produce these mosaic images, which have 25" pixels, 
is described in the text, (a) The Ha map, a combination of 89 individual observations, is 
presented with linear scaling, (b) The [O III] map includes 62 observations and is scaled 
linearly, (c) The linearly-scaled [S II] map is the result of 54 integrations, (d) A false-color 
combination of the three mosaics uses [S II] in red, Ha in green, and [O III] in blue. Blends 
of Ha and [S II] alone are yellow, [O III] and [S II] together appear as magenta, and Ha 
and [O III] combine to make cyan. The individual frames were median-filtered and scaled 
logarithmically before being combined into the false-color image. 

Fig. 2. — These mosaic images of portions of the Cygnus Loop are 95' across and have 6" 
pixels, (a) The northeast in Ha, (b) [O III], (c) [S II], and (d) a false-color combination as 
in Fig. (d). 

Fig. 3. — The northwest of the Cygnus Loop, as in Fig. 0. 
Fig. 4. — The southeast of the Cygnus Loop, as in Fig. |. 
Fig. 5. — The southwest of the Cygnus Loop, as in Fig. ||. 
Fig. 6. — The Cygnus Loop breakout, as in Fig. ^. 

Fig. 7. — Many of the features described in the text are noted on this Ha image. 

Fig. 8. — The true blast wave edge, detected in Balmer-dominated filaments, is nearly 
symmetrical. It is best fit by a circle centered on a = 20 h 51 m 21 s ,<5 = 31°01'37" with 
a radius of 1?4. Here the model circle is drawn on the ROSAT High Resolution Imager 
map of the Cygnus Loop in soft X-rays (Levenson et al. 1997). Although the X-ray 
observations are incomplete, having gaps on the interior near a = 20 h 49 m , 5 = 30°00' and 
a = 20 h 48 m , 5 = 31°30', they illustrate the blast wave clearly. 

Fig. 9. — This cartoon illustrates the orientation of a cloud along the line of sight near the 
roots of the "carrot," at a = 20 h 49 m ,<5 = 31°00'. As the blast wave expands radially (from 
ti through ti to £3), it interacts with the cloud from west to east in the plane of the sky, and 
the observer views its apparent progress along the arrow. 

Fig. 10. — Views of the interstellar medium before the supernova occurred. The three 
cartoons represent slices in the plane of the sky on the far side, through the center, and on 
the near side of the cavity, with north at the top and east toward the left. Features that 
appear in all three slices extend along the line of sight. A large wall of gas associated with 
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the Galactic plane and the cloud that causes radiative emission at the northeast are toward 
the rear. The atomic shell is prominent in the slice through the center of the cavity. Clouds 
at the northwest and around the edge of the breakout are on the near side. Scoville et al. 
(1977) observed the molecular gas near the Cygnus Loop. The extended wall of material at 
the west is detected in IRAS sky survey and in HI observations (DeNoyer 1975). The other 
clumps interact with the blast wave and are evident in the optical observations presented 
here and in the X-ray observations of Levenson et al. (1997). 
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